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COMPKESSmUJTY AND HEATmG EFFECTS ON PBESSUBE LOSS AND COOLING OF A 

BAFFLED CYLINDER BARREL 

By Arthttb W. Goldstein and Herman H. Ellerbrock, Jr. 



SUMMARY 

Theoretical investigatiom have shovm that, because air is 
compressible, the pressure-drop reguirements for cooling an 
air-cooled engine wHl be much greater at high aUitvdes and 
high speeds than at sea level and low speeds. Tests were 
conducted by the NAOA to obtain some experimenial con- 
firmaiion of the effect of air compressiMlity on cooling and 
pressure loss of a baffled cylinder barrel and to evaluate various 
methods of analysis. The results reported in the present 
paper are regarded as preliminary to tests on single-cylinder 
and muUicylinder engines. Tests were conducted over a wide 
range of airflows and density altitudes. 

The results indicate that, for a given air weight flow, the 
reduced pressure drop based on average cooling-^ir density 
Appap/ps, which has been used to correlate heatr-iransfer cool- 
ing data, is not constant for different air densities (Ap, cooling- 
air pressure drop across the engine; pa^/pt) ratio of average 
cooling-air density to density at sea level). Engine-cooliTig 
variables should therefore not be plotted against pressure drop. 
From the present tests a correlating variable for heatriranefer 
data is shown to be the air weight flow; the reduced pressure 
drop is not suitable for this purpose. An analysis based on 
the assumption of uniform flow is shown to be satisfactory for 
estimating the effect of compressiMlity on data obtained in these 
tests. A simpler empirical method in which compressibility 
and heating effects can be estimated was found for correlating 
the test data on pressure loss. 

INTRODUCTION 

Some investigators have heretofore correlated cooling data 
for air-cooled engines with the cooling-air pressure drop^ 
(references 1 to 4). The effect of compressibility was taken 
into account by using the product of the pressure drop and 
the average of the air densities at the front and the rear of 
the cylinder as the correlating variable rather than the pres- 



sure drop alone. With the high rate of heat exchange and 
the high air velocity between the fins required for effective 
cooling at high altitudes, however, a large air-density change 
will result. This change in air density is attended by an 
increase in velocity, and an additional pressure loss will be 
incurred at the baffle exit where this momentum will be 
lost. Estimates of the increase in pressure loss caused by 
air compressibility in engine-cooling systems with baffled 
cylinders were made in references 5 to 9. The analyses of 
references 5 and 9 were based on the assumptions of one- 
dimensional gas dynamics, but no experimental data were 
available to support these analyses. 

An investigation was begun by the NACA to obtain 
experimental confirmation of the effect of compressibility 
and rate of heat transfer on pressure loss and cooling of a 
baffled section of a cylinder barrel and to evaluate the 
various methods of correlating these data. The tests 
covered a range of simulated density altitudes from 4000 
to 33,000 feet, of velocities between the fins from a Mach 
number of 0.05 to near sonic values, and of heat inputs 
from 0 to 500 Btu per hour per squaxe inch of cylinder- 
wall smface. This investigation was conducted at Langley 
Memorial Aeronautical Laboratory, Langley Field, Va., 
during 1942. 

Acknowledgment is made to Mr. Frank E. Marble of the 
Supercharger and Airflow Research Division, Aircraft En- 
gine Eesearch Laboratory of the NACA for his suggestion 
that the Prandtl-Glauert compressibiiLty factor be used to 
correlate the results of the present tests. 

ANALYSIS 

ONE-DEVIENSIONAL GAS DYNAAHCS 

The analysis of the flow aroimd a baffled cylinder is based 
on the assumptions of one-dimensional gas dynamics. The 
equations for the pressure and density changes through a 
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baffled-cylinder system (fig. 1) are developed by the follow- 
ing analysis, which is similar to that of reference 5. 

1. The heat transfer between stations 1 and 2 and between 
stations 2 and 4 is estimated. 

2. No loss in total pressiu^e is assumed from station 1 to 
station 2. 

3. The pressure and the density at station 2 are then cal- 
culated from the heat-transfer estimate, the stagnation 
pressure, and the mass flow. 




FiauBB L— Finned cylinder, baffles, taflplece, and measuring-station locations used in tests. 

4. The analysis of the flow between stations 2 and 3 is 
based on the assumption of uniform flow of a compressible 
fluid with friction in a straight duct. 

5. The loss at the baffle exit is computed from the mo- 
mentum equations for uniform flow at stations 3 and 4. 

Change in gas state between station 1 and station 2. — If 
the heat picked up by the air between station 1 and station 
2 and the weight of air flowing through the baffle are known, 



where 



(1) 



rate of heat transfer to air between stations 1 and 2, 
Btu per second (The method of obtaining Hi is 
given in appendix A.) 



W weight of air flowing through baffle, pounds per 
second 

Cj, specific heat of air at constant pressure, Btu per 

poimd per ° F 
2^2. t stagnation temperature at station 2, ^ F absolute 
Ti,( stagnation temperature at station 1, ° F absolute 

A complete list of the symbols used is presented in appendix 
B. 

Stagnation or total temperatme and pressure as used in 
this report indicate gas properties that would obtain if the 
kinetic enei^ of the moving gas were isentropically con- 
verted into enthalpy. The temperature Ti.i may be con- 
sidered equal to the static temperattu-e at station 1 because 
the velocity at station 1 is negligible. If Tj,;, W, and Hi 
are known, the total temperature at station 2 may bo 
obtained from equation (1). 

The assmnption is made that there is no loss in total 
pressure from. station 1 to station 2 and that the stagnation 
state of the air at station 2 is known. Then, from the 
relations 

W 

*^=P2T^, (2) 



and 



gA2 



(3) 



and by means of the relations for isentropic change, the true- 
stream density and pressiure may be eliminated and the 
stagnation density and pressure inserted to give 



0 



7— 



IB 



(4) 



where 

acceleration of gravity, feet per second* 
cross-sectional area of free-flow space at station 2, 
square feet 

density of air at station 2, slugs per cubic foot (based 
on P2 and Tj) 

velocity of air in fln passage at station 2, feet per 

second 
Mach number at station 2 
ratio of specific heats for air (1.3947) 
static pressure of air at station 2, pounds per square 
foot 

stagnation air pressure at station 2, pounds per square 
foot 

stagnation air density at station 2 (computed from 

t and <), slugs per cubic foot 
true air-stream temperature of coohng air at station 2, 
absolute 



g 



P2 

M2 

T 

T2 
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FiouBG 2.— Helatton between ataguatlon state and static state of air as determined by mass 

flow. 

When {p%V^^l'p%,xP2,i is known, can be determined from 
equation (4); the ratios p^/pa, t flJid ^^a/^^a, < can be determined 
from the isentropic-change relations. A chart based on 
equation (4) and the isentropic relations is given in figure 2. 

Change in gas state between station 2 and station 3. — The 
pressure drop through the baffle from station 2 to station 3 
is given by the following equation, which has been modified 
from the corresponding equation in reference 5 by inserting 
the value of Qjy,f in terms of Ci?,/,*: 

Vmi^%lA ( i+Bl)+2(^-l) (5) 

(Z2 2 V P3/ \P3 / 



where 



dynamic pressure at station 2, Pal^a^^, potmds per 
square foot 

friction-drag coefficient between stations 2 and 3 

based on 22 (jJv, f^^^ 

friction-drag coefficient between station 2 and sta- 
tion 3 based on average of and 



friction drag from station 2 to station 3, pounds 

dynamic pressure at station 3, psT^s^^, poxmds per 
square foot 

static pressure at station 3, poimds per square foot 
density of air at station 3, slugs per cubic foot 
velocity of air at station 3, feet per second 

If the value for the local friction-drag coefficient Od is 
assiuned constant for all elements of the path, D/ would 
actually be obtained by a process of integration. The value 
of (7d should therefore be calculated from Df by means 
of some integrated mean value for JipFj. The value calcu- 
lated from the arithmetical average of 22 and gs, however, is 
used as an approximation. If this approximation is good, 
tiius calculated should be independent of inlet-density 
variations for fixed values of Eejniolds number and equal to 
the value that would be obtained with an incompressible 
fluid. By means of the continuity equation, the relationship 
between Cx)./ and Oj>./,< can be estabUshed as 



(6) 



The density ratio required for the solution of equation (5) 
can be calculated from the following equation, which is 
modified from the corresponding equation in reference 5 by 
inserting the value for (To,/ in terms of Od, f,t to give 



P3 



T-1 



'Hi) 



where 



S3 rate of heat transfer to air between stations 2 and 3, 
Btu per second (The method of determining 
is given in appendix A.) 

A convenient method of determining p2lpz uses the variables 



T' 



p'=Pt/pz 

p'-l 



v [%^* (i+p')+(p'-i)]+ V (1-"") 



(8) 
(9) 

(10) 

(11) 



When use is made of these definitions, equation (7) becomes 

Ma^=M^+rM,^ (12) 

Because Ma and M^,^ are functions only of and Ci,, i, for 
each value of Oj},f,i the variable Ma may be plotted against 
Mi,^ and p' may be plotted against by use of the defini- 
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tions for Ma and Af^*. These plots are shown in %ure 3. The 
value Ma must, however, be the ordinate satisfying the linear 
equation (12) with M^^ as the y intercept, M^,^ as the abscissa, 
and T' as the slope. The intersection of this line with the 
curve of Ma^ against Mj,^ for the given ODj.t determines 
Afa^. From the cmrve of Ma against p' for the given Gd j.i the 
value of can be obtained. 

Change in gas state between station 3 and station 4. — The 
momentum equation characterizing the pressure loss between 
station 3 and station 4 is derived for the baffle with a tail- 
piece. This loss as given ia references 5 and 9 was for a 
baffle without a tailpiece, where station 4 was in a section of 
very large area with very low velocity. For the present 
case the momentum equation is 

W 

2?3^ sin ez-p,A^+ fpdS^-Dz= -j- (V^- sin ^3) (13) 



where 

^3 angle between radii of cylinder to cylinder rear and to 

station 3, degrees (See fig. 1.) 
P4 static pressure at station 4, pounds per square foot 
^3 cross-sectional area at station 3, square feet 
Ai cross-sectional area of baffle exit at station 4, square 

feet 

j> static pressm:e, pounds per square foot 

dSx projection of any element of cylinder-wall surface and 
of curved part of baffle siuiace bacS: of station 3 on 
plane of Aa, square feet 

2>3 component of drag force normal to plane of A4 and ef- 
fective between ^stations 3 and 4, poimds 

V4 velocity of air at station 4, feet per second 
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The integral of pdSx is taken over the entire surface bounding 
the fluid between stations 3 and 4 except at flow cross sections 
of stations 3 and 4. 
It can be shown that 



fdSx=Ai—Ai sin flj 



(14) 



The pressure gradients will be proportional to the kinetic 
energy at station 3. Therefore, if 



/ (p-y3)<^» 



(15) 



if the coefficient for friction drag in the baffle exit is 

_ A 



0,. 



(16) 



and if these symbols and the continuity equation axe used, 
the momentum equation may be expressed as 



(17) 



In this equation the term 2{AzlA^ sia ^3 indicates that the 
momentum loss caused by the fact that the stream at the 
rear of the cylinder is not directed straight back cannot be 
neglected. This direction of the stream, however, tends to 
increase the value of Oz and the pressure recovery. The 
recovery coefficient 03, therefore, in some measure compen- 

sates for the fact that ^35*^90^ When-^ is small, all the 

terms become small and consequently pz is approximately 
equal to p^. 

For the present report the losses across the cylmder will be 
represented by equations (5) and (17). In the actual prac- 
tice of predicting the pressure losses, once the coefficients 
of equation (17) are evaluated, the density change may be 
evaluated from equation (17) and the energy equation. 

APPLICATION OF PBANDTL-GLAUERT COMPRESSIBHITY FACTOR TO 
FLOW ACROSS A BAFFLED CYLINDER 

The foregoing theory for flow of a compressible fluid is 
applicable only to the case of uniform velocity or one-dimen- 
sional fl ow. Fo r two-dimensional flow, the Prandtl-Glauert 
factor -^JI—Mq is used to compute the effect of compressi- 
bility. The Prandtl-Glauert factor is strictly applicable in 
flow conditions quite difi'erent from those existing around 
the baffled cylinder now being considered; therefore, this 
name should not be used for this application of the factor 



but will be used for convenience. The proper application 
of this factor is to static-pressure variations in the flow field 
of a body causing small perturbation velocities in an infinite 
uniform flow field with a frictionless, compressible fluid and 
with no heat transfer. In that case the difference in static 
pressure at any point in the flow field from what the static 
pressm^e would be with an incompressi ble fluid can be com- 
puted by the Prandtl-Glauert factor -yjl—Mf? in the equa- 
tion (reference 10) 



<^P.Jf^ 



or 



Op M — " 



(18) 



(19) 



where 

pressure of fluid that is characteristic of flow, pounds 
per square foot 
2o djmamic pressmre of fluid that is characteristic of flow 

^^PoT^o^^, pounds per square foot 

Mo Mach nimiber that is characteristic of flow 

Vi static pressure at same point as p with incompressible 

fluid, pounds per square foot 
Po density of fluid that is characteristic of flow, slugs per 

cubic foot 

velocity of fluid that is characteristic of flow, feet per 
second 

Ov,M pressure-loss coefficient with compressible flow 
Ojt^i pressure-loss coefficient with incompressible flow 

The factor -yjl—M,? was used to reduce the total-pressure- 
loss data to values that might be expected without compres- 
sibility effects. The application of the factor to the present 
data is to be regarded solely as an empirical method of cor- 
relation. No rational basis for the use of this factor is 
presented herein. 

The pressiu*e drop used in place of J?— i?o hi equation (18) 
for the baffled cylinder is Pi^t—PA.u The method of calculat- 
ing go and Mo was empirically determined by finding the 
pressm*e and the temperature that would give the best 
correlation of the data with different densities and various 
rates of heat transfer. The factor Op,i was then computed. 



(20) 



The total-pressme loss to be expected with an incompressi- 
ble fluid Lpi is then 



(21) 



843110—50 13 
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For the same mass flow (pV) but with standard density 
p,, Op^i will remain constant because it depends on only the 
Reynolds nxmaber. Then 

Apt^,=ApiPQ/p, (22) 

where 

Api^ loss in total presure under standard density condi- 
tions, pounds per square foot 

p, standard density (at 29.92 in. Hg and 60*" F), slugs 
per cubic foot 

From equations (20), (21), and (22) 

Api, ,=Aj?, ^ VI -Mo* (23) 
P$ 

where 

^Pt=Pi.t—P4,t 

THE BFFECT OF COMPRESSIBILTrT ON HBAT TKANSFER 

The effect of compressibility on the heat-transfer coeffi- 
cients is estimated by the effect calcidated for a flat plate 
by an equation given in reference 11. If the temperature 
of the plate is assumed to be 300® F (an average fin and 
cylinder temperature to be expected with a rear spark-plug 
temperature of 450° F), the Reynolds number is assumed to 
be 3000, and the free-stream air temperature is assumed to 
be —67® F (the temperature at high altitudes where com- 
pressibility effects in air-cooled engines may become critical), 
the effect of the free-stream Mach number M on the local 
heat-transfer coefficient hx is given by 

where Ax.i is the local heat-transfer coefficient that would 
exist with an incompressible fluid. The effect will be slightly 
less than shown in the foregoing equation for two reasons: 
(1) the value Ax.* is proportional to the local skin-friction 
coefficient, which in reference 12 is shown to decrease shghtly 
with increase of Mach number; (2) the highest possible value 
of AP= 1 will not exist over the entire cylinder. In practical 
cases M} will not be imity at any point aroimd the cylinder. 
Consequently, the effect of compressibiUty on heat-transfer 
coefficients for the cylinder can be expected to be negUgible. 

APPARATUS AND TESTS 

TEST SETUP 

The copper-plated steel barrel section used in the tests 
was 1% inches long with a 5Ji-inch bore. The fins on the 
cylinder were 0.50 inch wide, 0.036 inch thick, and spaced 
0.105 inch. Inside the cylinder was a grate to aid the 
pickup of heat and to reduce the temperature variation 
aroimd the inside of the cgrlinder. A metal baffle with a 
6-inch tailpiece was fitted around the cylinder (fig. 1). 
The tailpiece allowed the cross flow at the back of the 



cylinder to diminish sufficiently to permit more reliable 
pressure reading. The cross-sectional area of the exit of 
the baffie was 1.6 times the free-flow area between the fins. 
The unit was placed in an asbestos-lined metal box (fig. 4) 
and sealed at all edges with furnace cement to prevent air 
leakage. 

The source of heat was an oil biurner with a capacity of 
1 or 2 gallons of oil per hour, depending on the burner nozzle 
used. A firebrick furnace provided the space for the com- 
bustion to be completed before the hot gases came into 
contact with the cylinder grate. An auxiliary blower was 
needed with a large nozzle to supply the necessary air for 
combustion. 

The flow of cooling air was created by two compressors 
used as vacuum pumps and operated in series. Each pump 




FiQUBE 4.~Test unit showing test cylinder, baffles, taUpIoce, Insulated motal box, prossuro 
tabes, and thennocoaples. 



was driven by an 85-horsepower engine. In order to pre- 
vent suiting and to smooth out the flow of cooling air, largo 
tanks were placed upstream and downstream of the test 
cylinder. Bleed valves were placed in front of and between 
the pumps to provide fine control of the air flow. A tank 
with thin-plate orifices in each end was placed upstream of 
the test cylinder to measure the quantity of cooling air. A 
throttle placed between the orifice tank and the upstream 
surge tank was used to control the pressure of the cooling 
air in front of the test cylinder. A diagram of the apparatus 
is shown in figure 5. 

INSTRUMENTATION 

Surface temperatures of the cylinder were obtained at 24 
points (fig. 6) by means of insulated 28-gage chromel-almnol 
thermocouples spot-welded to the steel. The cold junctions 
of the thermocouples were inserted in a sealed wooden box. 
The temperature in the box was obtained with an alcohol-in- 
glass thermometer; the thermocouple potential was measured 
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with, a potentiometer. Thermocouples were also used to 
measure the orifice temperatures and the cooling-air inlet 
and outlet temperatures. The cooling-air outlet tempera- 
ture was measured downstream of the baffle tailpiece in an 
expanded section; it was therefore unnecessaiy to correct 
the readings for air velocity. The outlet duct was lagged 
to prevent heat loss. The accuracy of the temperature 
measurements was within ±1° F. 

Total and static pressures were obtained in the tailpiece 
and at two stations on the cylinder; a total pressure was 
obtained immediately in front of the cylinder (fig. 6). The 



Cooling- air 
intake 

Throtth valve 



Combustion -air 
i n fake ^^^^ 



■■o 



Furnace 



AuxlUory blower Oil burner 




Bleed vo/ves 



Pump 

Cooling -air exfyaust^^ 
Figure 5.— Diagram of test setup. 



total-pressure tubes of 0.030-inch-diameter steel tubing with 
a 0,006-inch-diameter hole in the side of each tube were 
inserted vertically around the cgrlindei-. The static pressures 
around the cylinder were measured by means of "vertical 
tubes inserted through holes tapped into the fins (fig. 7). 
Station 2 was at the baffle constriction and station 3' was 
somewhat behind the baffle-expansion point. Because of 
the location of the tube, the static-pressure readings for 
station 3' were not used in the calculations. Also, in view 
of the small distance between 3' and 3, the reading 'Pz',t was 
used for p3,i. Conventional-type pressure tubes were not 
used between the fins because it was thought they would 
block too much of the channel. A rake of conventional 
total-pressure tubes was used in survejdng the total pressure 
in the tailpiece, and a wall tap was used to measure the 



static pressure (fig. 1). The pressm-es wei-e read on vertical 
water or mercuiy manometers, depending on the range of 
pressures being measured. 

The very small pressure drops across the thin-plate 
orifices were measured with a micromanometer. 



Bo/f/e Ink 



Sfafion 2 




ThermocoupiG 



-Sfafion 3' 



o Tofal'prassure fube 
o 5f of ic' pressure fube 
riQUEE 0.— Location of thcnnocouples and pressure tnbcs on test cylinder. 



Tofal- pressure fube 



J] 



X////////. 



' Sialic- pressure fube 



Total- pressure top - 



\ ////////A , 



Static-pressure fop - 



7 



VtouSB 7.— fikatch sbowlng method of Installing pressoro tabes between fim. 




Test 
c\/linder 
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TESTS 

Teste were conducted over a range of air flows at several 
NACA standard density altitudes from 4000 to 33,000 feet 
simulated at tlie front of the cylinder. At each altitude, 
data were obtained with and without heat transfer for Mach 
numbers ranging from low values to the highest values 
obtainable with the apparatus. The highest Mach numbers 
occurred at the baffle esit, but no accurate data were ob- 
tained at that station. In the following table are given the 
ranges of Mach numbers at the baffle entrance and at 
the baffle exit ikfo; the ikfo values were computed on the as- 
sumption of no friction up to that point and are therefore 
lower than the maximum Mach nxunber of the flow system. 



Cylinder condition 


Density 
altitude 








f 4,000 
16.000 
23,000 
33,000 
4,000 
HOOO 
24,000 
32.000 


0.06 to 0.68 
.12 to .70 
A7to J58 
.18 to .60 
.16 to .73 
.16 to .67 
.21 to .69 
.20 to .66 


0.05 to 0.68 
.12 to .70 
.17 to .68 
.18 to .60 
.18 to .77 
.17 to .69 
.22 to .74 
.31 to .69 


With heat transfer 





Pressures and temperatures were recorded only after the 
test cylinder had practically reached a state of thermal 
equilibrium. The maxunmn allowable rate of change of 
temperature was about 3° F in 5 minutes for the maximum 
temperature before readings were taken. The cooling-air 
outlet temperature was read before and after reading the 
cylinder temperatures and an average of the two air tem- 
peratures was \ised. The two readings in no case differed 
by more than 2 percent of the temperature rise. Tempera- 
tures were read to within 0.4*^ F and pressures to within 
0.01 inch of mercury or water. The accuracy of the pressure 
readings was limited by the unsteadiness of the engines, 
which caused pressure fluctuations of as much as 0.1 inch 
of mercury. 

Summaries of the reduced test data and the derived 
quantities without and with heat transfer are given in tables 
I and n, respectively. Cylinder temperature-distribution 
data are available upon request to the NACA. 

COMPUTATIONS 
From the teste the following quantities were determined : 

^3, jA.2y jA^f -At, Alff and St 



where 



A, 



area of cylinder wall at base of fins, square feet 
average wall temperature, °F absolute (arithmetic 
average of measured temperatures of cylinder at 
base of fins) 

average temperature of cooling surface (fin and 

barrel surfaces), °F absolute 
total heat-transfer surface of cylinder, square feet 
From these data, the fictitious exit density pex and the 
coefficiente A, U, Cuj^ii and 03— C$ were determined as 
follows: 

Compulation of h. — The average surface heat-transfer 
coefficient h was calculated from the equation 



St 



h= 



The average temperature T/^av of the fin and barrel surface 
was obtained from the measured surface temperatures by 
averaging the temperatures with weighting factors propor- 
tional to the area elemente in which each thermocouple 
was located. 

Computation of U. — The average wall heat-transfer coef- 
ficient U was calculated from the following equation 



(24) 



Computation of fictitious exit density p^. — The stagnation 
pressure at station 3', the mass flow, the exit stagnation 
temperature, and the cross section at station 3 were used to 
compute the static pressure at station 3 by means of figure 2. 
This pressure and the temperature are used to compute 
the density p^x by means of the general gas law. 

Computation of Cd,/,<. — ^By the method given in appendix A, 

— TfT^ may be computed from the cooling-surface tom- 

perature distribution and the over-all rate of heat transfer. 
On the assumption that P2,t=Pi.tf Pt.t and (pVy/v^,t Pa.t 
were calculated and, from figure 2, pa, Pi, and ^2 wore 
found. From ^3., and Tz,t=Ti^tj the values of pz, qz, and 
Pz were calculated with the aid of figure 2. A value for 
OD,f,i was then computed from equation (5). 

Computation of exit coefflcients. — The exit coefficiente C3 
and Oz could not be separately determined and the value of 
Oz—Oz was consequently calcidated from jg, ^3, and pz 
(obtained in the calculation of Cz>./,<) and equation (17). 
From and ^4,^, figiure 2 was used to find the needed 
values for and p^. 
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Computation of Oj}j,ihj estimation of 03— C3. — ^The value 
of Gd,/,< was also calculated by the method of reference 6; 
that is, by estimating the exit-loss coefl5cient and assuming 
Ps=pi instead of using the measured value of 2?3,/- Equa- 
tion (17) was then used to calculate pz. Also from previous 
computations pa, fe, T\ and ilfa*, were known. The 
chart (fig. 3) was recalculated for constant values of 
(Pfl— 3?3)/ff3 instead of for constant values of Onj^i- The 
new chart was used to find the value of p' from (pa— i'8)/22> 
, iWj^, T^; equation (5) was then used to compute Ci>.r.<. 

Computation of Reynolds number. — ^Reynolds numbers B 
were obtained from the formula 

where 

d hydraulic diameter of fin passage, feet 

/i absolute viscosity of air based on average of average 
cooling-surface temperature and average cooling-air 
temperature, slugs per second per foot 
Computation of Op.i. — ^The pres§ure4oss coefficients Cp.i 

were calculated from the data by means of equation (20), 

which was used in the equivalent form 

n (1— p^/M.r ) Vut Po_ fTZTicn 

2 LPo.i Po. J 



From the data it was foxmd that po.t should be calculated from 
Put and Ti^t and that yo./ should equal pi,t for the best cor- 
relation. In the foregoing equation po,i is the stagnation 

density. The factor was then calculated and, from 

the value of this factor, was determined from 

' P0.f 

figure 2. 

RESULTS AND DISCUSSION 

EVALUATION OF COEFFiaBNTS 

Correlation of heat-transfer coefflcients h and V. — ^The 
efTect of altitude on the heat-transfer coefficients h andU" is 
shown in figures 8 and 9. The curve for h, the surface heat- 
transfer coefficient, plotted against weight flow of air at iihe 
baffle entrance shows no effect due to altitude inasmuch as 
the cmves for 4000 and 14,000 feet bracket the spread of the 
data. The curve for Z7, the wall heat-transfer coefficient, 
shows the same characteristics as the curve for h. These 
figures show that the Mach number had no effect on the heat 
transfer for the range covered. An effect of heat load was 
noted in the data but is not shown in figures 8 or 9. The 
heat-transfer coefficients h and U increase slightly with 
heat-transfer load. 
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Correlation of pressure-loss function (A^Par)/K(pTO*- — ^The 
pressure-loss function {^ppav)l}i{pV)^ is a dimensionless form 
of the conventional pressure-loss f miction Appajp,. The 
density is the average of densities upstream and down- 
stream of the cylinder, and p, is the air density xmder stand- 
ard sea-level conditions. When {Appap)/)i(pV)^ is plotted 
against weight flow pVg as in figure 10, the conclusions drawn 
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(a) Without heat transfer. 

(b) With heat transfer. 

FiQxmE 10 —Variation of pressuro-loss functbn with cooling-air weight flow. 

from the plot of the pressure-loss function may be applied 
to the function Ajppa,/p,. The dimensionless form has the 
advantage of showing compressibility effects more clearly. 

The plot of the pressure-loss function £^ppavl%{p^^ is 
shown iu figure 10 for nms with and without heat transfer. 
These curves show disagreement for different altitudes and 
for a variation in heat transfer. Theorise shown on the high 
weight-flow region of each curve indicates an additional 
pressure loss caused by compressibility effects at high Mach 
numbers. The high Mach numbers occur at lower weight 
flows as the density altitude iucreases. Any correlation 
method must make it possible to correct for these pressure 
rises and the separation of the curves for different altitudes 
and rates of heat transfer. 



Correlation of pressure-loss function (typpt^lVi{pYf, — 
Some investigators have proposed correlating the pressure loss 
by means of the fxmction Ajppcr, where p^x is the air density 
at the exit of the cylinder. This method was unsuccessfully 
tried with the data of the present tests. In most cases where 
this method would be apphed, the baflBie exit would have no 
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(a) Without heat transfer. 

(b) With heat transfer. 

Figure 11.— Variation of pressure-loM function with coolIng*Qlr weight flow. 

tailpiece and consequently no pressure recovery at the baffle 
exit. The exit density would then be computed from the 
static pressure at the baffle exit and the stagnation tempera- 
ture at the tailpiece. The density thus computed is lower 
than any density existing in any part of the test rig used for 
the present tests. The pressure-loss fimction {Appta^DiipV)'^ 
was computed by means of the density p„ and is plotted in 
figure 11 for tests with and without heat transfer. Because 
Pex decreases as compared with pat for increasing pressure 
losses, the correlation of the pressure-loss function 
(A^>p«)/}4(py)^ is better than that of (Ai?p«,)/K(pTO^. Al- 
though the correlation of the runs with heat transfer (fig. 
11(b)) is good, the runs without heat transfer (fig. 11(a)) 
show compressibility effects iu that the data for the different 
density altitudes form separate curves. 
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Correlation of drag coefficients and exit-recovery coeffl- 
cients. — ^The drag coefficients Cdja computed from the test 
data by the method previously given are plotted in figme 12 
against the weight flow per unit area. The utter lack of 
correlation between the results for the tests without heat 
transfer (fig. 12(a)) and those with heat transfer (fig. 12(b)) 
or between results at different altitudes is evident. The 
same result is to be observed in figure 13 in which the exit- 
recovery coefficient Oz^Cz has been plotted against the 
cooling-air weight flow. 

In order to find the source of error in these results, the 
assumptions involved in this method were examined. "With 
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(a) Wfthout heat tamafer. 

(b) With beat transfer. 

FiODBE 12.— Variation of baffleMshamicI drag coefHdent with cooHng-afr weight flow. Drag 
coeffidont calcoJated from test data. 



regard to the assiunption that the total pressure remains 
unchanged between stations 1 and 2, the data indicated that 
the pressure drop is so small as to be neghgible; at high flows 
the pressme drop amounts to about 1 percent of the total- 
pressure loss. 

The mass flow was calculated from the total and static 
pressures at stations 2 and 3' to test the validity of the 
assumption that the flow is uniform across each section (one- 
dimensional flow). The calculation was made from the 
factor {pVyiftPh which was obtained from the pressure 
ratio 2?/^?; and figme 2. The ratios of the weight flows to 
the weight flows calculated from measurements at the 
orifice are plotted for stations 2 and 3' in figm-es 14 and 15, 
respectively, against the weight flows obtained from the 
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FiGUBE 18.— Variation of exit coefficient with oooling-air weight flow. 
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orifice measurements. At station 2, figm*e 14 shows that 
the ratio is near enough unity for appHcation of the assump- 
tion of imiform flow. At station 3', however, figure 15 
indicates that the measurements are not accurate at low 
flows. In the calculation of any coeffident hased on measure- 
ments at one point at station 3', correlation will probably 
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FiQxmB 14.— Batlo of mass flow measorod al station 2 to mass flow at station 2 obtatned from 
orlflce mcdsaremoits. Tests made withoat heat transfer. 

not be obtained. It is partly for this reason that CTdj-,!, 
a splotted in figure 12, shows such lack of correlation. 

In order to eliminate the use of these inaccurate measure- 
ments, Czjj-.f was recalculated by the use of an assumed 
exit loss instead of the measmred loss. The value of 03 
compensates to some extent for the fact that Qzj^WP. If 

03+2^ sin ^3 is assmned constant or various values of 

^a, then for ^3=90° separation from the cylinder rear can be 
expected and no pressure recovery due to pressure gradients 
along the wall will result. In that case 03=0. 
Then, 

a,+2^sm53=2^ 



or 



08=2-^ (1— sin^a) 



(26) 



If it is further assumed that C^=0 andj( pa=P4 and if 
q^—qj^iA^iAzf- is substituted in equation (17), then 



This equation appUes with an incompressible fluid flowing 
through the sudden expansion of a straight duct. 

The method of calculating Cz>. r. < when this loss is used has 
been previously given. The results are plotted in figure 16. 
Much better correlation results tidan wap pbtained in figure 
12, especially at low flows. The data for each altitude form 
a smooth curve. The curves separate, however, as soon as 
the Mach number, which increases with altitude for a flxed 
Reynolds nxmiber, becomes appreciable, indicating that the 
effect of compressibihty has been overestimated by taking 
too high a loss at the exit; that is, 03 is too low. 

In order to determine how much these irregularities in 
CTij.i and 03— C3 affect theover-aU total-pressure drop, the 
values of Odj.i and 03— determined from the tests were 
used to calculate the loss in total pressure to be expected 
with an incompressible fluid by assmning that the density 
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FiOTTRK 15.— Ratio of mass flow measured at station 3' to mass flow at station 3' obtained from 
orifice measurements. Tests made withoat heat transfer. 

remains unchanged in all equations for pressure loss. If aU 
losses are added and if p= constant =P2, 

0,.,^^^,= l+i7z>./.i+a3-a3-2 ^ sin53+ (0 (26) 

The over-all pressure-loss coefficient (7p,i is plotted against 
the weight flow per unit area in figure 17, which shows that 
the extreme scatter of the data points in figm^es 12 and 13 
has been sufficiently reduced to allow fairing of curves 
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(a) Wlthoat beat tiaiufor. 

(b) With heat traiufer. 

FiousE 10.— Sffeot of oooUng-alr iirel^t on IHctloB-diag ooefSdant calculated by i 
of assumed estt loss. 



through the data points. The effect of compressibility may 
be seen from the fact that the curves show systematic 
differences for the different altitudes although the ends of 
the curves do not have the sharp curvature that is evident 
in the plots of {^fPa^l%{pVY. These differences are much 
smaller than the entire compressibility effect, as may be 
seen by comparison with figure 10, which is already some- 
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(a) Without beet transfer. 

(b) With heat transfer. 

Piaims 17.— Effect of oooUng-aIr weight fow on pressoie-Ioss coefficient calculated iiom test 
data by analytical method. 



what corrected for compressibility effects. This result 
would seem to indicate that, for correction of the compressi- 
bility effect, the division of the pressure loss between the 
baffle channel and the baflSe exit is unimportant because 
the pressure was no doubt in error. In order to investigate 
this hypothesis, the over-all pressure-loss coefficient Cp.< 
was computed from the estimated value for 03— Cs (equation 
(25)) and the corresponding values for Goj^t (%. 16). 
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The results of these calculations are plotted in figure 18, 
which shows that the value for cannot be estimated in an 
arbitrary fashion. 

If the values of Ci>./,< and (h—Cs from equations (5) and 
(17) are inserted in equation (26) in order to determine how 
the pressure-loss coefficient C^.f depends on the data, there 
is obtained 

^Kt)"('-?>(0 

If the density ratios are estimated from the energy equa- 
tions for very low Mach mmibers (M=0), then 



Pt—Ps 



2 



2» 1^(1+2") 



2T'_ 
T' 



The value of T' (equation (8)) is always less than 0.10 and 
the data indicate that the term involving is never greater 
than 0.01. At low Mach numbers, therefore, large varia- 
tions in estimates of pz from its true value will not aflFect 
(7p.<. For high Mach numbers, however, equation (27) 
indicates that pz will definitely aflFect the result for (7p,f. 
Figure 15 shows that the measurements of pz are inaccurate 
at low Mach niunbers but much better at high Machnimibers ; 
therefore, at high Mach nmnbers when the data for pz are 
used, fau* correlation of Op,< may be expected. At low 
Mach numbers, the inaccurate values for jpz do not aflFect the 
result; correlation may consequently also be expected in this 
range. With inaccurate values for (Tdj-.* and Oz—Oz at low 
Mach nmnbers, correlation for Cj,,i may be expected; 
whereas, in the high range, correlation may be expected only 
when accurate values for Cij,/.i and Oz—Cz are obtained. 

Figure 12(a) shows that no correlation of Cd,/.< occurs at 
low Mach numbers and fair correlation occurs at high Mach 
numbers. No accurate correlation in any range is obtained 
for tests with heat transfer (fig. 12(b)) but, if scatter is taken 
about a mean curve for the high Mach numbers, the varia- 
tions are small compared with the value of Op^i and therefore 
do not show up so prominently. The same general features 
of the Ci>, Afplots (fig. 12) may be observed in theoa— C3 plots 
(fig. 13) and the same comments apply. From the C7),/.< 
data (fig. 16) it may be presumed that the estimate of 03 was 
incorrect and resulted in the calculation of faulty values for 



fz' Comparison of Cp^i of figure 17 with that of figure 18 
indicates that, in the range of low Mach numbers for all 
altitudes, the C7p,< curves are the same with the measured and 
the calcxilated values of pz) which confirms the deduction 
that, for low values of Mach number, Cp^is independent of 
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(a) Withoat heat transfer. 

(b) With heat tranafer. 

FiQUBE 18.— Effect of cooUng-elr weight flow on pressure-loss coolllolont colcubtod by luo 
of assumed exit loss. 

fz' In fact, the lowest-altitude curve is approximately the 
sajne for both methods of computing Cp,i except at the very 
highest Mach numbers. The lack of correlation of the data 
of figure 18 at high Mach nxmabers confirms the deduction 
that the estimated values of pz are incorrect and influence 
the values for Cp^i and that the measured values for high 
Mach nmnbers used in calculating (7p,< (fig. 17) are accurate. 
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The fact that the pressure measurements in the rear of the 
cylinder were veiy inaccurate at the low cooling-air weight 
flows indicates a lack of uniformity in the flow, which con- 
tradicts the hypothesis upon which the flow analysis was 
built. A second inaccuracy- is seen in the hypothesis by 
which Cd f iis calculatcfl from Cd,/- The method of calcula- 

increases linearly along the channel; J^fe+ffa), therefore, 
represents a good mean value for q. This assumption is 



tions is based on the assimiption that g_ 
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FiauRE lO.^EfTcct of Reynolds number on pressure-loss coeiBdent computed by empirical 

method. 



approximately true for low Mach nimabers but not for high 
Mach numbers where a much more rapid rate of increase is 
to be expected. Because ^(22+23) is higher than the mean 
effective value of 2 in the channel, Cz?,/,< must be lower than 
the coefficient that may be expected with a truly incompres- 
sible fluid. The nonuniformity of flow in the rear (station 3) 
\vill also modify the momentum equations for the baffle 
channel and the baffle exit. For these reasons values of 
Ct?, Ai and 03—63 obtained by the assumption of uniform flow 
and by means of equations (5) and (17) will probably not 
check estimates from measurements of the velocity and 
pressure distribution along the channel and baffle walls- 
Further experiments are required to determine conclusively 
the effectiveness of this method. 

Correlation of the coefficient (7p i calculated by empirical 
method. — plot of the pressure-loss coefficient Oj,^i against 
the Reynolds number is shp\vn in figure 19. Although the 
correlation shown in this plot is satisfactory, some featwes 
show that improvement is still desirable. The points for 



the lowest density altitude without heat transfer are, in 
general, lower than the average of the points. The tendency 
for this plot to faU off at high Reynolds nimibers indicates 
that the correction for compressibihty effects is probably too 
large. This result might possibly be due to the fact that 
the correction factor apphed is too large for the pressure 
losses occurring before station 3. This effect does not occur 
in the data for other density altitude^s, possibly because the 
data for the lowest densities may be inaccurate in the high- 
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FiQUBK 20.— Variation of corrcctod prossaro loss with cooling-air wofght flow. 

flow region as a result of the imsteadiness of the driving 
engines. 

In order to determine the magnitude of these deviations 
on the actual pressures, the corrected pressure drops were 
plotted in figure 20 in the customaiy manner. The dis- 
crepancies noted in figure 19 show little effect on the pressure- 
loss plot. The shape of the plots of figure 19 indicates that 
a linear curve of pressure-loss data is not to be expected. 
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GENERAL DISCUSSION 

Comparison of tlie correlation methods. — ^For any set of 
pressure-loss data, the separation of the curves for different 
density altitudes is an effect of compressibility that has not 
been corrected. The effect of compressibiHty is accentuated 
at the highest Aveight flow for each density-altitude curve. 
The effectiveness of each method of pressuxe-drop correlation 
is characterized by the magnitude of the spread of the curves 
at the four weight flows that are the maximums for the dif- 
ferent altitude curves. At the end points of the curves for 
4000-, 14,000-, 24,000-, and 32,000-feet density altitudes, 
there are, respectively, two, four, six, and eight points on 
all the curves for comparison. The average deviation from 
the mean of all the curves at these points is given iu the 
followiDg tables ia percentage of the mean value of the para- 
meter at that poiat: 





Mean deviation, percent 


Weight flow 








Ob/(sec) 




Method 




(«!")) 












Analytical 


14,4 


13 


7.3 


7.8 




13 


3.9 


4.7 


20.0 


11 


2.8 


6.0 


4ao 


13 


4.9 


2.1 



Reynolds 
number, R 


Mean devia- 
tion, empir- 
ical method 
(percent) 


ai4X10» 


4.8 


.21 


2.4 


.30 


2.1 


.40 


L5 



These tables show that the use of the variable {Appav)l%(pVy 
is least effective, that the use of the variable (App<x)/K(pT^)^ 
and the analytical method is better and about equally 
effective, and that the empirical method using the Prandtl- 
Glauert factor is most effective and of satisfactory accuracy. 

The involved procedure needed for use of the analytical 
method is a definite drawback. Purthermore, certain iofor- 
mation is needed for pressure-loss predictions that is not 
necessary when the empirical method is used; the relative 
amount of heat picked up by the air in front of the baiffle 
entrance must be obtained and pressures at the baffle-channel 
expansion point are required in order to evaluate separately 
the coefficients Cj>./.< and at—Cz in the pressure-loss predic- 
tions. The testing technique for determining reliable values 
of pz is yet to be developed. Computations made for a 
modem engine at an altitude of 40,000 feet show that very 
highMach numbers occur and that separate values of 
and (h—Oz are needed. The method is somewhat simplified, 
however, when no baffle tailpiece is used, which elimiaates 
pressure recovery at the baffle exit. In that case P3=Pii and 
pressure measurements are therefore unnecessary in the 
baffle exit. 



The same drawback applies to the use of the pressure-loss 
fimction (Appej)/K(pTO^; that is, a measurement of the 
pressure in the baffle exit is required to determine p« unless 
a baffle is used that gives complete loss of the Mnetic energy 
at the baffle exit. This method of correcting the pressure 
loss for compressibility effects has a rational basis. If the 
entire loss is assumed to occur at the baffle axit, if no pressure 
recovery exists in the baffle tailpiece, and if tiie air temper- 
ature in the baffle exit is assmned equal to that behind the 
cylinder, {App^)/%{pV)^ should very nearly account for all 
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FiGUBE 21.— Prediction of coollng-elr pressore loss at altitude by xmot ^(pV)t ^ -^i— jVfo>, 

assumed a function of Reynolds number, and by use of y^^y^t » assumed a function of 
mass flow. 



compressibility effects by equation (17). If not all the loss 
in energy occurs at the baffle exit and if some pressure 
recovery is present there, this method should give good 
results if most of the pressure loss occurs at the baffle exit 
and if the rest varies in much the same manner as the 
baffle-exit loss. 

The empirical method of computation of the pressure 
loss has the advantage of being quite simple. There is 
doubt as to the types of flow apparatus to which the empirical 
method can be apphed because no rational basis for it 
exists at present. Test data from engines, engine cylinders, 
and radiators are required before an estimate of the reliability 
of this method can be made. -The same confirmation of 
the reHabUity of the analytical method is needed, however, 
because there is some doubt that it correctly represented 
the actual flow conditions in even this comparatively simple 
type of flow path. 

In order to illustrate the application of the empirical 
method, an example is computed and plotted in figure 21. 
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For comparison, the pressure loss is also computed by means 
of the function (^pPav)l%{pV)^* The barrel temperature is 
assumed constant at 350° F; the rate of heat transfer is 
assumed constant and such as to require an air weight flow 
at sea level of 18.6 poxmds per second per square foot of 
free-flow area of the barrel-baffle channel. The atmospheric 
pressures and temperatiures are modified for isentropic 
compression corresponding to a flight speed of 200 miles per 
hour. The subsequent steps are: 

1. Compute cooling-temperature differential Tb— Ti,,. 

2. Compute the heat-transfer coefficient U (equation 
(24)). 

3. Determine pVg (fig. 9) ; compute the exit temperature. 

4. Compute the Reynolds number, and determine Op < 
(fig. 19). 

5. Compute the abscissa of figure 2 and read the compres- 
sibility-correction factor (po/po,f) -yjl-^M^. 

6. Compute the pressure loss from the data of steps 3, 4, 
and 6. 

When the pressure loss is computed by means of the 
function {^fPa^l%{pVYf the first three steps are the same 
as those previously given. The weight flow is used to 
determine (Ai?PaF)/K(pV)^ by means of the plot of lowest 
density altitude (fig. 10 (b)). An estimate of Aj) is made for 
determining p^,; the resulting value of A^? is used to refine 
the values for pat and Aj?. For less extreme cooling condi- 
tions, the difference between the two methods is not so great. 

Effect of viscosity, — ^For the general correlation of pressure- 
loss and heat-transfer data, dimensionless parameters shoifld 
be used. The variable for example, should be plotted iu 

terms of the Stanton number against the Rejniolds 

number. In turbulent fiow the Stanton number varies 
approximately as the viscosity to the 1/5 power. For prac- 
tical appHcation of heat-transfer data to altitude-cooling 
predictions, therefore, the small variation in the viscosity 
caused by the variations in temperature of the cooling-air 
surface film will affect the heat transfer only sUghtly. Good 
correlation of heat-transfer data of a given engine or cyl- 
inder can be expected, therefore, if the weight fiow is used 
as the correlation variable. The same reasoning can be 
applied to pressure-loss data. ' Because the present pressure- 
loss data were obtained with and without heat transfer, 
considerable difference should be expected in the film 
viscosity. Consequently, the Reynolds number was used 
to correlate the pressure-loss data; the viscosity of the 
surface film gave better correlation than the cooling-air 
viscosity in the computation of the Reynolds number. The 
pressure-loss coefficients were plotted against the weight 
flow because use of the Reynolds number effected no increase 
in correlation or significant change in the plot. 

Velocity distribution. — Some improvement might be made 
in the analysis of flow in radiators, for which the assumption 



of imiform flow at the radiator-tube exit is satisfactory, if 
allowance were made in the momentum equation for the 
nonmuformity of flow to be expected with a tm-bulent veloc- 
ity distribution in the tube. The computed kinetic-enei^ 
loss at the exit would thus be increased. 

SUMMARY OF RESULTS 

Based on an analysis of pressure-drop requirements and 
on experiments to determine the effect of air compressibility 
on cooling and pressure loss of a baffled cylinder barrel, the 
following results were obtained: 

1. The pressmre loss from the front of the cylinder to the 
baffle entrance was very small, as might be expected for any 
baffled engiue cylinder. 

2. The method of analyzing the flow processes around a 
baffled cylinder based on the assumption of imiform flow 
corrected for most of the effect of compressibility. 

3. The assmnption of uniform flow used in the analytical 
method of predicting pressure losses was not verified by 
computations from pressure measurements between the fins 
at low flows. 

4. Prediction of the pressure losses across a cyliader at 
high Mach numbers usiug the analytical method based on 
uniform flow requires knowledge of the exact values of the 
friction coefficient between the fins and the coefficient at the 
exit of the cylinder. At low Mach numbers the division of 
the pressure coefficients has httle effect on predicted pressure 
loss. 

5. The use of a fictitious exit density p^ in correcting 
pressure-loss data gives as accurate correlation as the analyt- 
ical method based on the assmnption of uniform flow. 

6. An empirical method Jias been found which gives satis- 
factory correlation (mean deviation of pressure-function 
curves of about 3 percent) of the test data on pressure loss of 
the present cylinder and which permits estimates to be made 
of compressibflity and heating effects more simply than either 
the analytical method or the method using the baffle-exit 
density in baffle-flow systems with a tailpiece. 

CONCLUSION 

If engine-cooling equations are based on the cooling-air 
weight flow instead of on the pressure loss, good correlation 
of data and prediction of performance may be expected 
because test data and theoretical estimates showed no appre- 
ciable effect of air compressibility on heat-transfer coeffi- 
cients in the ordinary range of engine operating conditions. 



Atbcraft Engine Research Labobatobt, 
National Advisory Committbb for Abronautics, 
Cleveland, Ohio, Jvly 1, 1944, 
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ESTEVfATION OF AIR-TEMPERATURE DISTRIBUTION AROUND A BAFFLED CYLINDER BARREL 



In order to determine the pressure loss over the cylinder, 
it is necessary to know the rate of heat transfer to the cooling 
air between stations 1 and 2 Hi and the rate of heat transfer 
to the cooling air in the baffled section S2. The following 
analysis presents one method of estimating Hi and H2. 

For an element of heat-transmitting surface dS, the amount 
of heat transferred per unit time is 

dH=h(Tf-Tt)dS (Al) 

where 

H rate of heat transfer from front of cylinder to point being 

considered, Btu/sec 
hg local surface heat-transfer coefficient, Btu/(sec) (^F) (sq ft) 
Tf cooling-surface temperature average at any local flow 

section, °F absolute 
Tt local stagnation-air temperatin^e, ^F absolute 

The viscous shearing stress slows down the air in the 
immediate vicinity of the wall and increases the temperatxure 
there until the temperature upon which the heat transfer 
depends, is more properly the stagnation temperature than 
the static temperature (reference 11), If Wl is the local 
weight flow between the fins (lb/sec), the enei^ equation is 

H=c^W,,(T,-T,,,) (A2) 

In the unbaffled section of the cylinder, T7i increases from 
a value of zero at the front stagnation poiat to the value W 
at the baffle entry. The term Tt is eliminated from equa- 
tions (A2) and (Al) and the quantities 

and 

^sT^— 2\ ,=cooling-temperatm'e differential, °F 
are used to obtain 

dH=^ G(cpWLt'-n)d{S/ST) 

where 

St entire surface for heat transfer, sq ft 
S surface for heat transfer from front of cylinder to point 
considered, sq ft 
The solution of this differential equation is 

£r=exp [-J^^^^(7rf(S/Sr)]j Constant-!- 

rsisr ' r rsisr n ) 

\ CcpWx^i exp Cd{SIST)\ diS/Sr) } (A3) 

Because S/St=^0 and ^=0 at the front of the cylinder, the 
constant of integration is zero. At the baffle entrance, 
202 



S/St is equal to a, the ratio of the unbaffled cooling surface 
of the cylinder to the total cooling surface, and H—Hx, the 
rate of heat exchange in the unbaffled portion of the cylinder. 
At the baffle exit, S'/St'=1, and n=Hi+Hi, the rate of heat 
exchange over the entire cylinder. When these boundary 
conditions are inserted in equation (A3), 

fli=exp j^Cd(S/ST)j 

jy<h>WLi exp[J^^^'ca(iS/Sr)] diS/Sr) (A4) 

Hx+H,=exp ^-j^CdiS/Srij 

£cCpWz.i exp^j^^^Od{S/Sr)J diS/Sr) (A5) 

The rest of this appendix derives expressions for Hi and 
H1+H2J which are more convenient than equations (A4) and 
(A5) for computation from the test data. The quantity C 
varies little over the cylinder except at the very front. The 

factor J^^'^^^^Vd the mtegral of CWjei^^^'^^^^^^ d(SISr) 
are only slightly affected by this variation. Consequently, 
the quantity (7 is considered a constant in aU the integrations. 
The values found for Hi and H1+H2 from the energy 
equations 

H^=Wc,{T,,,-T^,,) 

and 

are substituted in equations (A4) and (A5) to obtain 

Ti.t- Tr,,= Ce-^'j^ ^ ie<^(sisr)d (S/St) (AO) 

£'^feW;d(S/Sr)] (A7) 

Ce-^P ^ ic^W^ri diS/Sr) (A8) 

If an estimate is made of the -nreight-flow variation Wi./W 
and measurements of t are obtained, C can be determined from 
the over-all temperature rise 14,1— Ti,, and equation (A7) 
and its value used in equation (A6) to find Ta,,— Ti.t. 
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In order to estimate the weight-flow variation in the nn- 
baffled portion of the cyUnder, the flow in the front is com- 
puted \)j means of the equation for potential flow over a 
circular cylinder. The velocity over such a cylinder is 

where Vq ia the velocity of the undisturbed flow, which 
reaches a maximum at S/St=' 1/2, In the case of the baffled 
cylinder the maximum value is at S/St=oc. The velocity 
function is therefore modified to take this factor into account 

^=2^0 sin (irSl2ccST) 

The constriction at the baffle entrance {S/St=cc) of the 
flow area to I/ti its value for an unbaffled cylinder increases 
the maximum velocity by the factor n, but the velocity in 
front of the cylinder is unchanged. The velocity distribution 
must therefore be multipHed by an even function of S/St, 
which increases from a value of 1 at S/St=0 to a value of n at 

S/St= a. The assumed function is j^l — ^^^^ cos ('n-fi^/2aiS'r) J 



which gives an approximate weight-flow variation of 
-^=sin (irSl2aST) sm J (S/St) 



(A9) 



From the baffle dimensions, n=5. For the unbaffled portion 
of the cylinder, a linear cooUng-surface temperature dis- 
tribution is assumed 



(AlO) 



where li is the value of t at S/St=Oj and is the value of 
i at SISt=cc, Equations (A9) and (AlO) are substituted in 
equation (A6) to obtain 

(.T,, ,-T,, ,)/f,=/i(C)+(^)/,(C) (All) 

where 



n—1 



{20|+,-[2<7l+,(0.+(l^ 



An approximation must now be found for the integral 
involved in equation (A8) . In the baffled section, Wl/W= 1 



and t varies very Uttlo; an extreme case shows 9-percent aver- 
age deviation from the mean. The mibaffled cylinder rear 
shows a much larger deviation but the decrease in the effec- 
tive cooling-air velocity reduces the variation of the function 
(Wi,/W)t, The value of T7i/T7is therefore taken to be unity 
and an average value used for t so that 



j: 



^ te^3/Sr)d(^S/ST) =^ (e^-6c«)f3 (A12) 



where ^3 is the average cooliog-surface temperature of the 
cylinder from the baffle entrance to the rear minus the inlet- 
air temperature. This result is substituted in equation (A8) 
to obtain 



If 



(A13) 



equation (All) can be plotted in the form (la.*— ^'i.O/^a as a 
function of (e^^^"**^— 1) and iti—h)/^' This plot is shown 
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in figure 22. Equation (Al3), however, gives a solution for 
the variable {T2^t^Ti,t)/ta as a linear function of e^^^'^^—l 
with a y intercept of (Ti,t^Ti,t)fta and a slope of 
[iz—(T4,t'-'Ti,t)]/ta' From data on the temperature of the 
cylinder, the straight line (equation (Al3)) can be located on 
figure 22 and the rutersection point can be obtained with the 
ciu^e of equation (All) for the proper value of (ti—i2)/ta; 
thus, (T^^t^Ti^^Jta can be determined. 



APPENDIX B 



SYMBOLS 



A cross-sectional area at a station indicated by a sub- 
script (except with '^subscript b), sq ft 
At area of cylinder wall at base of fins, sq ft 
Cp specific heat of air at constant pressure, Btu/(lb) CF) 

C3 coeflScient for friction drag in baffle exit 

CjDj- coefficient for skin-friction drag from station 2 to 

station 3 based on dynamic pressure at station 2 
Cj>^.i coefficient for drag from station 2 to station 3 based 

on dynamic pressure averaged in baffle channel 
Cp,i pressure-loss coefficient with an incompressible fluid 
<7p,jf pressure-loss coefficient with a compressible fluid 
g acceleration of gravity, ft/sec^ 

kx local surface heat-transfer coefficient, Btu/(sec) (**F) 
(sq ft) 

h average surface heat-transfer coefficient, based on 

inlet-air temperature, Btu/(sec)C'F)(sq ft) 
H rate of heat transfer from cylinder area between front 

and any point considered, Btu/sec 
Hi rate of heat transfer to air from front of cylinder to 

baffle entrance, Btu/sec 
H2 rate of heat transfer from baffle entrance to baffle 

exit (station 2 to station 4), Btu/sec 
M Mach number at a station indicated by subscript 

(except Ma, M^) 

" ■\l 7P'[%i (1+P') + (P'-1)]+V 

p static or stagnation indicated by subscript t) pres- 
sure at any point in fluid indicated by a subscript, 
Ib/sq ft 

Pt static pressure at any point in fluid with incompres- 
sible flow, Ib/sq ft 

2 dynamic pressure at any point in fluid indicated by 
subscript, Ib/sq ft 

R Reynolds number 

S surface for heat transfer from front of cylinder to 

point considered, sq ft 
St total heat-transfer surface of cylinder, sq ft 
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Sx projected area of cylinder-wall surface and curved 
part of baffle surface back of station 3 on plane 
of A4J sq ft 

tz average value of t between stations 2 and 4, ° F 
T true air-stream or stagnation (indicated by subscript i) 
temperature at any point in fluid indicated by 
subscript (except Tt,, Tf, and T/.a,), ° F absolute 
average temperature of cyhnder wall at base of fins, 
F absolute 

Tf cooling-surface temperature averaged at any local 

flow section, ^ F absolute 
Tf^av average temperature of cooling surface of cylinder, 

° F absolute 
T s Bj/cpTTTa 

U wall heat-transfer coefficient of cylinder, Btu/(sec) 
(sq ft) r F) 

V velocity of air at any station indicated by subscript, 
ft/sec 

W weight of air flowing through baffle, lb/sec 

Wl local air weight flow between fins, lb/sec 

a ratio of unbaffled cooling surface of cylinder to total 

cooling surface 
7 ratio of specific heats for air (1.3947) 
^3 angle between radii of cyHnder to cylinder rear and 

to station 3, deg (See fig. 1.) 
II absolute viscosity, slugs/(sec)(ft) 
p local true or stagnation (indicated by subscript 0 

density at any point indicated by subscript (except 

p„ p„, and p«), slugs/cu ft 
p, standard density of air at 29.92 in. Hg and 60** F, 

slugs/cu ft 

Pap average of densities at stations 1 and 4, slugs/cu ft 
fictitious exit density 
s P2/P3 

Ap total-pressure drop from front to rear of cylinder, 
Ib/sq ft or in. water, (pi,(--^>4.i) 

Lpi loss in total pressure from front to rear of cylinder 
with an incompressible fluid, Ib/sq ft or in. water 

A^f,, loss in total pressure from front to rear of cylinder 
with an incompressible fluid under standard den- 
sity conditions, Ib/sq ft or in. water 

Subscripts apphcable to A, Pj T, T, p, M, and g: 

1, 2, 3, 3', 4 stations indicated in figure 1 

i stagnation condition of gas (absence of subscript t 
indicates true-stream condition) 

0 condition characteristic of entire flow 
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TABLE I.— SUMMARY OF DATA FOR TESTS WITHOUT HEAT TRANSFER 
[The symbols used are defined in appendix B] 
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36.89 


20.27 
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.4926 
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1.021 
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1.0620 
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18.88 


.637 


.622 


.0296 


-.2611 


.6216 


.6737 


.2260 


&46 


18.02 


.680 


.649 


.1479 


1405 


.6138 


.6839 


.2680 


11.21 


15.91 


.811 


.716 


.4478 


-1839 


.6383 


.6712 


,2770 


14.87 


16.36 


.846 


,712 


.4661 


.1923 


.6294 


.6551 


.3060 


17.82 


13.80 


.941 


-742 


.4362 


.1724 


.6384 


.6554 


.8340 


21,25 


12.76 


.965 


.712 


.4987 


.1604 


.6639 


.6494 


.3363 


21.37 


19.39 


.695 


.682 


.4314 


.1743 


,6451 


.6776 


.1228 


2.83 


19.26 


.842 


.823 


.8333 


.1769 


.7448 


.8132 


.1461 


4-85 


18.81 


.777 


.747 


.3984 


.1606 


.6843 


.7356 


,1808 


6.© 


18.22 


.781 


-736 


.4266 


.1963 


.8719 


.7222 


.2100 


8.86 


17.80 


.761 


.693 


.4822 


,1888 


.6088 


.6781 


.2360 


ia39 


16.21 


.868 


.762 


.4456 


,1961 


.6617 


.n79 


,2913 


16.10 


14.35 


.998 


.981 


.1846 


.3430 


1.0606 


.9803 


,0807 


1.84 


13.99 


1.007 


-974 


.2883 


.2688 


.8827 


.9688 


.1286 


4.21 


13.67 


.898 


.857 


-2679 


.1416 


.7869 


.8458 


.1492 


6,08 


12.88 


.876 


.811 


.3570 


,1689 


.7141 


.7848 


.1827 


7,20 


12.13 


.890 


.791 


.3817 


.1383 


.6788 


.7606 


.2104 


9.23 


14.27 


1.010 


.995 


.0188 


-1703 


L0537 


.9738 


.0947 


2,60 


14.15 


.887 


.884 


,1877 


.2186 


.9331 


.8803 


.1077 


2.83 


14.01 


-875 


.860 


.2572 


,2690 


.9040 


.8304 


.1180 


3,34 


laTD 


-863 


.824 


.4816 


.1624 


-6031 


.8037 


.1366 


4,30 


13.16 


.877 


-828 


,3039 


,1894 


.7877 


.8046 


.1621 


6.02 


12.73 


.867 


.797 


.3685 


.2094 


.7531 


-7658 


.1813 


7.22 


1L83 


.892 


.784 


.3948 


.1771 


.6846 


.7379 


.2157 


9.56 


11.76 


-901 


.772 


.3862 


.1772 


.6942 


.7137 


.2253 


ia33 


9.62 


.988 


.774 


.4061 


.1521 


.8482 


.6826 


.2631 


1Z33 


10-81 


L089 


L072 


.0767 


.3262 


L1527 


L0690 


.0625 


1.20 


10.12 


.961 


.931 


.3279 


.3708 


.9461 


.9263 


.0743 


L48 


9.99 


.983 


.948 


.3390 


-3886 


.9197 


.9391 


.0880 


2.02 


7.26 


1.027 


.890 


.2167 


-1640 


.8095 


.7961 


.1683 


6.61 


7.96 


1.017 


.902 


.2168 


.1961 


.8617 


.8363 


.1639 


5.76 


8.81 


.987 


.912 


.2193 


.2231 


.8760 


.8874 


.1300 


4.30 


9.64 


L034 


-080 


.2995 


-3438 


.9163 


.9619 


-1061 


3,09 



Density altitude (ft) 



4,000.< 



15,000., 



23,000., 



33,000., 



abs.) 



28.10 
28.12 
28.19 
28.20 
28.34 
28,20 
28,11 
28,14 
28.20 
28.24 
28.16 
28.16 
28,14 
28.23 
28.22 
28.23 
2a 18 
28.09 
28.11 
28.16 
20.08 
19.93 
20.03 
19.92 
2a 12 
20. 11 
20. 18 
20,12 
20.04 
2a 09 
20. 16 
20.02 
20.18 
20.20 
20.20 
2a 15 
20. 24 
20. 18 
14.82 
15,00 
15.07 
16.01 
16.07 
14.86 
14.83 
14.91 
14.91 
14.88 
14.91 
14.90 
14.89 
14.85 
ia77 
ia73 
10.79 
ia76 

laTo 
ia68 

10.72 



Ti.t 
C»F abs.) 



642.8 
642.8 
643.5 
644.6 
644.5 
645.6 
646.5 
646.6 
646.6 
645,6 
640.7 
641.8 
54Z8 
543.5 
649.2 
647.9 
548.2 
547.0 
66L0 
667.0 
649.0 
650.0 
560.0 
65L0 
651.0 
644.0 
649.0 
649.0 
649.0 
649.0 
65a0 
655.0 
647.6 
547.8 
647.6 
648.6 
647.6 
647.8 
662.0 
639.8 
642.6 
546.6 

«a6 

649.0 
649.0 
650.0 
660.0 
660.0 
66L0 
543.0 
649.0 
547.0 
662.0 
662.0 
663.0 
66Z0 
553.0 
564.0 
663.0 



ab/&c) 

(sqft)) 



3.892 
4-030 
4.491 
5.464 
6.477 
7.621 
9.230 
11.487 
13.388 
15,791 
18.217 
19.810 
22.173 
24.673 
26.772 
29.736 
33,375 
38.74 
40.60 
39.75 
6.608 
7.627 
9.494 
ia411 
11.68 
15.79 
19.46 
22.22 
23,74 
28,31 
28.81 
29.08 
ia36 
12.34 
16.26 
17.72 
19.84 
24.60 
8.93 

laeo 

12. 41 

16.36 
17.76 
8.07 
9. 18 

iao9 

11.68 
13.86 
15.64 
18.09 
19.21 
2L52 
6.38 
6.39 
7.39 
14. 46 
13.23 
1L20 
9.05 



ao8 

.07 
.09 
.09 
.14 
.16 
.22 
.29 
.42 
.68 
.79 
.91 
1.12 
L38 
1.73 
2.22 
2,99 
4.61 
5.14 
4.97 
.20 
.26 
.36 
.42 
.62 
.91 
1.14 
L62 
2.26 
2,97 
4. 12 
4.44 
,34 
.69 
.84 
L16 
L40 
2.43 
.30 
.69 
.86 
1.33 
1.84 
.41 
.47 
.66 
.74 
L09 
1.37 
1.97 
Z27 
3.31 
.27 
.34 
.47 
2.11 
L71 
1. 15 
.78 



(in^g 
abs.) 



28.11 
2a 12 
2a 20 
2a 28 
2a 31 
2a 19 
2a 11 
2a 13 
2a 26 
2a 23 
2a 16 
2a 16 
2a 13 
2a 21 
2a 21 
2a 20 

2a 14 

2a 08 
2a 04 

2a 10 

20.08 
19.94 
2a 03 
19,90 



20.08 
2a 14 



2a 01 

2a 06 

2a 16 
2a 01 

20. 14 
20.22 

2a 17 

20. 14 



20.14 
14.82 
16.02 
16.06 
14.99 
16.07 
14.89 
14,82 
14,91 
14.91 
14-89 
14.89 
14.90 
14.88 
14.84 
ia76 
ia72 
ia76 
ia75 

ia67 
ia67 

10.70 



pi 
an.Hg 
abs,) 



2a 04 
2a 05 
2a 11 
2a 15 
2a 18 
2a 01 
27,81 
27.73 
27.65 
27.45 
27-18 
26.96 
2a 68 
2a 41 
25.89 
25.35 
24.38 
22.28 
2L71 
2L89 
19.84 
19.61 
19.66 
19.34 
ia92 
iaS7 
ia20 
17.61 
17-04 
ia27 
16.40 
14.92 
19.57 
19.44 
19.01 
ia65 
ia20 
ia97 
14.44 
14.20 
ia94 
13.30 
12.78 
14.38 
14- 19 
14, 15 
13.92 
la 51 
la 14 
12.47 
12.08 
1L21 
ia40 
ia23 

la 14 
a35 
a69 

9.24 
9.78 



(in. fig 
abs.) 



2a 06 
2a 21 

2a 16 

2a 25 
2a 30 
2a 17 
2a 09 
2a 08 
2a 10 
2a 06 
27.97 
27,90 
27.82 
27.81 
27-69 
27.62 
27,34 
27.02 
2a 97 
2a 98 
19.99 
19.96 
19.96 
19.86 
2a 00 
ia88 

20.73 
20.66 
19.46 
20.30 
19.23 
ia96 
20.06 
20.08 
2a 02 
19.84 
19.86 
19.47 
14.71 
14.80 
14.92 
14.72 
14-70 
14.80 
14.73 
14.83 
14.78 
14.68 
14.82 
14.44 
14.37 
14- 17 
ia70 

laoo 

ia62 

ia44 
ia38 
ia44 
ia67 



Pi' 

(in-Hg 
abs,) 



2a 07 
2a 21 
2a 18 
2a 22 
2a 22 
2a 08 
27.87 
27.78 
27.84 
27.40 
27.09 
2a 82 
2a 40 
2a 16 
26.64 
24.89 
23,81 
20.90 
2a 09 
20.27 
2a 34 
2a 14 
2a 06 
19.89 
19.89 
19.33 
19.66 
ia74 
17.21 
16.16 
14.71 
ia76 
20.04 
19.88 
19-42 
ia88 
ia46 
ia92 
15.45 
16. 10 
14.84 
14.08 
ia43 
15.39 
16.21 
15-18 
14.89 
14-42 
14.03 
ia26 
12.79 
11.34 
1Z02 
1L87 
U-77 
9.68 
ia20 
ia88 
11.44 
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[The symbols used are defined in appendix B] 



Density altitude at) 


(in^g 
abs.) 


Ti,i 
(*F abs.) 


(Wft)) 


»' 


r F abs.) 


(i|*j|g 
abs.) 


Pi 

abs.) 


(in.Hg 
abs.) 


(in^g 
abs.) 


V* 
(In.Hg 
abs.) 


(° F) 


(»F) 


4,000 

14,000 _- 

24,000 

• 

32,000 


27.80 
28.00 
28.04 
28.02 
27.99 
27.94 
28.15 
28.09 
28.01 
2a 47 
2a 47 
2a 47 
2a 48 
2a 42 
2a 46 

2a 61 

2a £3 

14.98 

14.97 

16.03 

16.06 

14.87 

14.84 " 

14.84 

14.92 

1L18 

11.18 

1L19 

1L18 

ILIS 

1L24 


66a6 
65L1 
662.1 
663.1 
663.6 
655.6 
56ai 
647.1 
553.8 
54a4 
641,0 
641.3 
542.1 
512.1 
642.1 
542.1 
64a2 
54a 6 
66L1 
651.6 
666.1 
667.6 
667.6 
666.6 
660.6 
637.6 
53a3 
53a6 
64a 4 
512.1 
542.8 


12.22 
14.97 
18.16 
2L62 
2a 30 
2a 53 
3&12 
30.62 
4a 83 

a65 
1L87 
14.54 
17.26 
ia68 
22.31 
2a 87 
2a 06 

8.20 
1L40 
16.41 
2L62 
ia09 
13.90 
16.86 
18.26 

8.56 

a76 
1L37 
12.22 
14.01 
14.37 


a54 

.70 
LOl 
1.39 
2.12 
2.76 
4.69 
6.76 
7.60 
.39 
.61 
.86 
1.22 
1.64 
2.18 
3.44 
4.61 
.62 
.89 
L43 
4.10 
.77 
L35 
1.80 
2.41 
.60 
.92 
1.28 
L66 
2.22 
2.66 


63a6 
628.6 
623.1 
62a6 
613.6 
6U.6 
697.1 
603.6 
694.6 
616.8 
606.2 
597.4 
693.2 
689.1 
686.6 
68a9 
681.9 
626.6 

6iai 
6oai 

602.1 
66a6 
632.6 
824.6 
622.1 
603.3 
604.3 

eoas 

698.9 
698.2 
608.2 


27.77 
27.97 
28.03 
2&01 
27.93 
27.93 
28.13 
28.02 
27.96 
2a 48 
2a 44 
2a 44 
2a 44 
2a 37 
2a 45 
2a 60 

15.00 
14.08 
14.99 

14.86 
14.82 
14-81 
14.89 
1L19 
1L16 
1L16 
1L16 
1L17 
1L25 


27.27 
27.24 
2a 92 
2a 43 
25.68 
24.86 
23.60 
22.13 
2L69 
2a 07 
ia78 
ia39 
18.94 
18.38 
17.82 
ia63 

iai6 

14.62 
14.09 
13.44 
12.26 
14.13 
13.49 
13.06 
12.66 

lase 
laii 

9. 76 
9.49 
8.90 
8.83 


27.60 
27.76 
27.74 
27.64 
27.42 
27.21 
27.06 
2a 88 
2a 79 
2a 38 
2a 30 
2a 24 

2a 12 
laoo 

19.34 
19.74 
ia47 
14.80 
14.67 
14.67 
14.35 
14.61 
14.47 
14. i3 
14.30 
ILOl 

iao5 
ia92 
ia88 
lasi 
laoo 


27.06 
2a 94 
2a 60 
26.85 
24.65 
23.46 
2a 85 
17.67 
16.86 
2a 61 
2a 28 

ia84 

1&32 
18.63 
17.42 
la 12 
14.18 
14.71 
14.06 
13.26 
a36 
14.16 
ia34 
12.62 
1L69 
12.10 
1L77 
11.38 

iao7 

ia25 
a96 


27.02 
2a 06 
2a 49 
25.90 
34,77 
23.62 
2L 13 
18.06 
ia46 
1&90 
ia62 

iao6 

18. 51 
17.70 

laoi 

16.14 
13.05 
14.24 
13.54 
1176 
8.75 
13.05 
12.84 
12.18 
IL 10 

iao8 
a67 

a 17 

a60 

7.67 
7.23 


Sia66 

302.17 
283,60 
268,39 
26a 09 
23a 69 
202,00 
20a30 

isa 18 

234.96 
217.77 
198.07 
185.02 
172,00 
16a 00 
151.63 
142. 12 
261. 16 
227.34 
2(M.24 
183.00 
308.47 
203.83 
m77 
23a 82 
208.02 
2ia69 
20L45 
197.11 
10L02 
189.42 


307.00 
305.65 
36L16 
34a 09 
32L76 
3ia63 
28a 25 
273.38 
24a 21 
27a 40 
266.00 
230.00 
227.09 
217. M 
21L 17 
197.20 
186.88 
289.00 
208,M 
24a 02 
227.26 
30a65 
317.38 
30ai3 
287.00 

mo2 

245.02 
237.40 
23a 01 
231.70 
228.90 


Density altitude (ft) 


U 

(Btu/(hr) 


h 

CBtn/(hr) 








App„ 






Cdjj 


(as^t- 
method) 


(empir- 
ical 
method) 


Boynolds 
number 

xio-» 


(In. HiO) 


4,000.. 

14,000- 

24.000- 

32,000.- 


a636 
.768 
.878 
L020 
L172 
1.271 
1.402 
L610 
L605 
.671 
.726 
.814 
.927 
1.016 
1.098 
L261 
L221 
.616 
.664 
.812 
L040 
. ffi^ 
.788 
.868 
.960 
.658 
.616 
.699 
.726 
.810 
.830 


aos6i 

.1052 
.1246 
.1482 
.1728 
.1897 
.2301 
.2620 
.2369 
.0763 
.0976 
.1127 
.1305 
.1467 
.1602 
.1883 
.1831 
.0680 
.0900 
.1126 
.1480 
.0834 
.1086 
.1176 
.1364 
.0741 
.0324 
.0944 
.0902 
.1122 
.1151 


385.0 
36a5 

34ao 

323.6 
302.6 

28ao 

268.6 
264.5 
226.6 

27ao 

254.0 
238.0 
222.8 
20a8 
22a2 
19a 4 
18a 8 
284.0 
262.3 
23a6 
22L1 
363.0 
308.6 
29L0 

27ao 

238,0 
242.9 
234.1 
232.6 
227.1 

22ao 


a8055 
.8372 
.8629 
.8624 
.9167 
.9371 
L0&48 
LW12 
L0687 
.6444 
.7087 
.7663 
.8416 
.8270 
.9977 
.9883 
1.0039 
.6926 
.6646 
.7474 
-7879 
.6912 
.7912 
.8454 
.8986 
.0290 
.653^ 
.©29 
.7097 
.7683 
.7788 


28a6 
224.3 
203.6 
18a8 
I6a7 
16L4 
127.8 
12a3 
112.2 
183.4 
166.2 
144.0 
12a2 

iia7 
iia8 
9a2 
9ai 

201.4 

mi 

15L1 

i2a5 

24L3 
197.2 
181.6 
166.9 
163.8 
162.6 
15L4 

14a 0 
14a 1 

138.4 


1.019 

.880 
.863 
.848 
.860 
.846 
.948 
1.077 
LOSS 
L1S8 
.986 
.841 
.857 
.866 
.878 
.931 
LOIS 
L 160 
LOll 

L168 
LlOO 
.977 
LOlO 
.988 
.932 
L098 
L080 
L213 
L182 
L240 


a 931 

.802 
.780 
.749 
.736 
.723 
.747 
.783 
.733 

LWl 
.918 
.773 
.768 
.764 
.732 
.739 
.730 

L061 
.898 
.781 
.785 
.966 
.844 
.837 
.779 
.838 
.967 
.916 

L026 
.961 

L049 


a 4672 
.4606 
.4118 
.4190 
.4522 
.4632 
.4386 
.3918 
.4399 
.1167 
.2830 
.4281 
.4544 
.4740 
.6463 
.4616 
.4816 
.2844 
.4033 
.4203 
.3936 
.4888 
.3868 
.4004 
.4217 
.3784 
.2402 
.2666 
.1614 
.2224 
.1354 


a 2347 
.1612 
.1834 
.1070 
.1469 
-1630 
.0904 
.0724 
.0633 
.0927 
.1844 
.1674 
.1609 
.1933 
.3178 
.1216 
.1039 
.2994 
.2862 
.1632 
.0710 
.3862 

.ino 

.1629 
.1176 
.1802 
.1825 
.1561 
.1260 
.0747 
.0869 


a6797 
.6039 
.6738 
.6002 
.6959 
.6920 
.6540 
.6828 
.6266 
.8792 
.8036 
.6316 
.6087 
.6216 
.6747 
.6622 
.6246 
.9172 
.7341 
.6361 
,5797 
.7098 
.6873 
.6647 
.5981 
.7040 
.8446 
.8006 
.8868 
.7515 
.8627 


a9299 
.7^ 
.7741 
.7381 
.7173 
.0980 
.6013 
.6862 
.6340 

L0368 
.9151 
.7847 
,7686 
.7638 
.7343 
.8877 
.6647 

L0698 
.8900 
.7372 
.6980 
.9621 
.8306 
.8166 
.7322 
.8302 
.9146 
.9161 
.9910 
.8890 
.9600 


a 1126 

.1386 
.1600 
.2024 
.2601 
.2818 
.3644 
.3006 
.4047 
.0827 
.1111 
,1442 
.1722 
.1983 
.2269 
.2747 
.2969 
.0779 
.1107 
.1606 
.2121 
.0923 
.1308 
.1601 
.1730 
.0842 
.0960 
.1120 
.1211 
.1392 
.1429 


6,432 

aow 

a 087 
13.331 
ia406 
23.808 
36.280 
42.074 
41.330 
2.008 
4.020 
a322 
8.806 
n.417 
14.200 
ia422 
2L037 
2.787 
4.600 
a 847 
12.041 
3.830 
a332 
&016 

aw7 

2.374 
3.620 
4.632 
6.780 
a 826 
7.762 



